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Positive-parity spin-i excitations of the nucleon are explored in lattice QCD. The variational 
method is used in this investigation and several correlation matrices are employed. As our focus is 
on the utility and methodology of the variational approach, we work in the quenched approximation 
to QCD. Various sweeps of Gaussian fermion-field smearing is applied at the source and at the sink 
of XiXi E^nd X1X2 correlation functions to obtain a large basis of operators. Using several different 
approaches for constructing basis interpolators, we demonstrate how improving the basis can split 
what otherwise might be interpreted as a single state into multiple eigenstates. Consistency of the 
extracted excited energy states are explored over various dimensions of the correlation matrices. The 
use of large correlation matrices is emphasized for the reliable extraction of the excited eigenstates 
of QCD. 
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I. INTRODUCTION 

One of the long-standing puzzles in hadron spec- 
troscopy has been the low mass of the first positive par- 
ity, = i"*", excitation of the nucleon, known as the 
Roper resonance A^*(1440 MeV). In constituent or va- 
lence quark models with harmonic oscillator potentials, 
the lowest-lying odd parity state naturally occurs be- 
low the N —h^ state (with principal quantum number 
iV = 2) 3, if whereas, in Nature the Roper resonance 

is almost 100 MeV below the N = i~(1535 MeV) state. 
Similar difficulties in the level orderings appear for the 
= 1"^ A* (1600) and i"^S*(1690) resonances, which 
have led to the speculation that the Roper resonance may 
be more appropriately viewed as a hybrid baryon state 
with explicitly excited glue field configurations 0, 01 or 
as a breathing mode of the ground state [1] or states 
which can be described in terms of meson-baryon dy- 
namics alone Q. 

Lattice QCD is very successful in computing many 
properties of hadrons from first principles. In particular, 
in hadron spectroscopy, the ground states of the hadron 
spectrum are now well understood Q. However, the ex- 
cited states still prove a significant challenge as they be- 
long to the sub-leading exponential of the two-point cor- 
relation function. Extracting excited states masses from 
these sub-leading exponents is difficult as the correlation 
functions decay quickly and the signal to noise ratio dete- 
riorates more rapidly. In baryon spectroscopy, there are 
many experimentally observed baryon resonances whose 
physical properties are poorly understood. Lattice QCD 
can provide theoretical input to solidify their identifica- 
tion. The first detailed analysis of the positive parity 
excitation of the nucleon was performed in Ref. |8[ us- 
ing Wilson fermions and an operator product expansion 
spectral ansatz. Since then several attempts have been 



made to address these issues in the lattice framework 
[20j , but in many cases no potential identification of the 
Roper state has been made. Recently, however, in the 
analysis of Refs. [3, [l^, [2l| a low-lying Roper state has 
been identified using Bayesian techniques. 

Another state-of-the-art approach in hadron spec- 
troscopy is the 'variational method' [H], [2^, which is 
based on a correlation matrix analysis. The identifi- 
cation of the Roper state with this method have been 
mixed. However, recently, in Ref. (2^ a low-lying Roper 
state has been identified with this approach employing a 
diverse range of smeared-smeared correlation functions. 
Our work there motivates us to investigate the several 
positive parity excited states using similar techniques but 
in a significantly more comprehensive manner. 

In this paper, the variational analysis used in Refs. [l^, 
[25! is explored more extensively. In particular, we con- 
sider 6x6 and 8x8 correlation matrices not only built 
from the xi interpolating field, but also incorporating an- 
other nucleon interpolator, X2, to extend the set of basis 
operators. 6x6 matrices are built up using XiAi ^^'^ 
X1X2 correlators, while the 8x8 matrices use the X1X2 
correlation functions, as discussed in the text. 

One of the goals of this paper is to investigate the 
high- lying positive parity spin-i excited states of the nu- 
cleon, such as Pii(1710 MeV) and Pii(2100 MeV), using 
larger correlation matrices. Incorporating the X2 interpo- 
lator with various numbers of smearing sweeps, enables 
us to explore more deeply the overlapping of different in- 
terpolators with the energy eigen-states. This will also 
prove the reliability of the discovery of the Roper res- 
onance 2j|. We demonstrate how improving the basis 
interpolating fields can split what otherwise might be in- 
terpreted as a single state into multiple eigenstates. 

This paper is arranged as follows: Section |TT] contains 
the general description of the extraction of masses with 
the introduction of different nucleon interpolating fields. 
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The lattice details are given in Section IIIIl the results 
are presented in Section ITVl and conclusions arc made in 
Section El 



II. MASS OF HADRONS 

The masses of hadrons arc extracted from two-point 
correlation functions using operators chosen to have over- 
lap with desired states. Let us consider a baryon state 
B of spin half, if we suppress Dirac indices a two point 
function can be written as, 



J2 \B,p',s){B,p',s\^I, 



(2) 



where B can include multi-particle states. The substitu- 
tion of Eq. (121) into Eq. Q yields. 



X B,p'.s 

(3) 



Gij(t,p) — ''^e ^P''^ {il\T{xiix)xj{0)}\^i)- (1) We can express the operator Xi (a;) as 



The operator Xi(0) creates states from the vacuum at 
space-time point and, following the evolution of the 
states in time t, the states are destroyed by the operator 
Xi{x) at point x, t. T stands for the time ordered product 
of operators. A complete set of momentum eigenstates 
provides. 



(4) 



where, P'' = P = {H, P) and P is the momentum op- 
erator whose eigenvalue is the total momentum of the 
system. Eq. ^ can now be written as. 



Gy(<,p) =E E e-■'P■^{^\e'''''x^i^)e-'''^\B,p',s){B,p\s\x,{m) 
= E E e-'^«*e---(^^-^^')(17|x.(0)|i?,p',s)(i?,p',s|x,(0)|17). 

X B^p',s 



(5) 



As we move from Minkowski space to Euclidean space. For the negative parity states one requires 
the time t ^ —it and the above equation then can be 

written as, , 



G^Jit,p)^ E e-^"'Sp,p,{^l\x^mB,p',s){B,p',s\xJmn) 

B,p',s 

"e~'^^*{n\x^mB,P,s){B,p,s\x,mn). 



B 



(6) 



The overlap of the interpolating fields x(0) and x(0) with 
positive and negative parity baryon states \B^) can be 
parametrized by a complex quantity called the coupling 
strength, \g± , which can be defined for positive parity 
states by 



{n\xmB-,p,s) = XB-- 



-l5UB-{P,s), (9) 



Pr- 



■UB-{P,s)"f5. (10) 



Here, Xb± and A^i are the couplings of the interpolat- 
ing functions at the sink and the source respectively and 
Mb± is the mass of the state B^. Eb± is the energy of 

the state i?^, where Eb± = Mg± +p^, and ub±{p,s) 

and ub± {Pi s) a-re the Dirac spinors. 



(r!|x(0)|i?+,p,s) = A 



B+^ 



B + 



E 



B + 



UB+iP,s), (7) 



{B+,p,s\xiOm=>^B+^ 



B+ 



E 



B+ 



ub+{p,s). (8) 



(11) 



Thus, Eq. (ini) contains a projection operator r± 



l'B± 



{p, s)u'^^{p, s), through which the contributions 



to the even and odd parity states from the correlation 
function can be obtained. For positive parity, this can 
be expressed as, 
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7.p + 



2^; 



B+ 



(12) 



-7.p + Afg- 



(13) 



and for the negative parity, 



By substituting the above Eqs. for the positive and neg- 
ative parity states in Eq. ([6]) we obtain, 

r ^ \^ \ \ -E„+tl-PB+ + Mb+ 



B+ 



2Ef 



(14) 



[2 a ] in the spatial dimensions is apphed through an iter- 
ative process. The smearing procedure is: 



V'^(a:,i) =^J^(a:^,a^')V^.-i(a^',i), 



(18) 



where, 



F{x,x') = (1 - a)Sx,x' + 7 y^[^A'(-^)'^^',^-i-A 



(19) 



where the parameter a = 0.7 is used in our calculation. 
After repeating the procedures N^-m times on a point 
source the resulting smeared fcrmion field is, 

V'JV.„(a:,<) =E^'^^'"(^'^')'^o(^''*)- (20) 



At momentum p = 0, Eb± = Mb± , a parity projection 
operator T± can be introduced, 



B. Variational Method 



r± = -(l±7o). (15) 

We can isolate the masses of the even and odd parity 
states by taking the trace of Q with the operators r+ and 
r_ . The positive parity state propagates through the 
(1,1) and (2,2) elements of the Dirac matrix, whereas, 
negative parity state propagates through the (3,3) and 
(4,4) elements. 

The correlation function for positive and negative par- 
ity states can then be written as, 

G±(t,0)-Tr,p[r±g,,(i,0)] 

The correlation function contains a superposition of 
states. The mass of the lowest state, Mq± can be ex- 
tracted at large t where the contributions from all other 
states are suppressed, 

G±(t,0)*^°°A±A±e-^V*. (17) 

A. Source Smearing 

The spatial source smearing [2^ technique is applied 
to increase the overlap of the interpolators with the 
lower lying states. We employ a fixed boundary con- 
dition in the time direction for the fermions by setting 
Ut{x,Nt) — OVx in the hopping terms of the fcrmion 
action with periodic boundary conditions imposed in the 
spatial directions. Gauge invariant Gaussian smearing 



The extraction of the ground state mass can be done 
straightforwardly. However access to the excited state 
masses requires additional effort. Here we consider the 
variational method [l^, [l^]. The variational method 
requires the cross correlation of operators so that the op- 
erator space can be diagonalised and the excited state 
masses extracted from the exponential nature of the di- 
agonalised basis. To access N states of the spectrum, 
one requires a minimum of N interpolators. With the 
assumption that only N states contribute significantly to 
Gij at time i, the parity projected two point correlation 
function matrix for p = Q can be written as, 

= (X]Trsp{r±(l]|x.(x)x,(0)|17)}) (21) 

X 

= ^ A?A^"e-"°*, (22) 

where Dirac indices are implicit. Here, A" and A" are 
the couplings of interpolators Xi Xj 3.t the sink and 
source respectively to eigenstates a = 0, • • • , (iV — 1). 
is the mass of the state a. The N interpolators have the 
same quantum numbers and provide an iV-dimensional 
basis upon which to describe the states. Using this basis 
we aim to construct TV independent interpolating source 
and sink fields which isolate N baryon states \Ba), i.e. 

N 

r = E<^- (23) 

i=l 
N 

'^"=E<^- (24) 

i=l 
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such that, 



III. SIMULATION DETAILS 



= (5„/j2"u(a,p,s), (25) 

{n\q}"\Bp,p,s) ^S^pz"u{a,p,s), (26) 

where z" and z" are the couphng strengths of 0" and 0" 
to the state 1-8^). Consider a real eigenvector which 
operates on the correlation matrix Gij (t) from the right, 
one can obtain [l^l , 



X?z"e- 



(27) 



For notational convenience, in the remainder of the dis- 
cussion the repeated indices i, j, fc are to be understood as 
being summed over, whereas, a, which stands for a par- 
ticular state, is not. Since the only t dependence comes 
from the exponential term, we can write a recurrence re- 
lation at time {t + At) as. 



(28) 



for sufficiently large t and t + At [25ll27l|. 

Multiplying the above equation by [Gij{t)]^^ from the 
left we get. 



[{G{t))-^G{t + At)]u°' = e 



(29) 



This is an eigenvalue equation for eigenvector u" with 



eigenvalue c" 



We can also solve the left eigen- 



value equation to recover the v°' eigenvector. 



<G,,(t + At) 



<G,,(t). 



Similarly, 



v'^iGit + At){G{t))-^] =e 



(30) 



(31) 



The vectors and i;f diagonalizc the correlation matrix 
at time t and t + At making the projected correlation 
matrix. 



<jfG,j{t)u'^. = J"'3^a-/3g-m„t^ 



. (32) 

The parity projected, eigenstate projected correlator, 
vfG^At)uf = Gt is then used to obtain masses of differ- 
ent states. We construct the effective mass 

G^(t,0) \ 



M^sit) = In 



G5(t + 1,0) 



(33) 



and apply standardised analysis techniques as described 
in Ref. [ia. 



We use an ensemble of 200 quenched configurations 
with a lattice volume of 16"^ x 32. Gauge field configu- 
rations are generated by using the doubly blocked Wil- 
son action in two coupling space (DBW2) [2^, An 
0(a)-improved fat link irrelevant clover (FLIC) fcrmion 
action [30[ is used to generate quark propagators. This 
action has excellent scaling properties and pro vides near 



continuum results at finite lattice spacing |3l| . The lat- 
tice spacing is a = 0.127 fm, as determined by the static 
quark potential, with the scale set with the Sommer scale. 
To = 0.49 fm [33 • In the irrelevant operators of the 
fermion action we apply four sweeps of stout-link smear- 
ing to the gauge links to reduce the coupling with the 
high frequency modes of the theory [Sj]- We use the 
same method as in Ref. |3J] to determine fixed boundary 



effects, and the effects are significant only after time slice 
25 in the present analysis. 

Eight different levels of gauge invariant Gaussian 
smearing [2^ (1, 3, 7, 12, 16, 26, 35, 48 sweeps cor- 
responding to rms radii, in lattice units, of 0.6897, 
1.0459, 1.5831, 2.0639, 2.3792, 3.0284, 3.5237, 4.1868) 
are applied at the source (at t = 4) and at the sink. This 
is to ensure a variety of overlaps of the interpolators 
with the lower-lying states. The analysis is performed 
on nine different quark masses providing pion masses of 
m^={0.797,0. 729,0.641,0.541,0.430,0.380,0.327,0. 295,0. 249 
} GeV. The error analysis is performed using the jack- 
knife method, where the x^/dof is obtained via a 
covariance matrix analysis method. 

The nucleon interpolators we consider in this analysis 
are 



{u"'{x)Gj5d''ix))u''{x), 



{u"'{x)Gd''ix))-f5u''{x). 



(34) 
(35) 



We use the Dirac representation of the gamma matrices 
in our analysis. 

Each of the 6x6 correlation matrices of XiXi corre- 
lators corresponds to a particular selection of 6 levels of 
smearing from the 8 that we have available. We consid- 
ered six different combinations of these to give six dif- 
ferent correlation matrices. Our selections are shown in 
Table |T1 We note that 48 sweeps tends to be noisy and 
therefore eliminate it from most of our considerations. 

For the 6x6 matrices of x 1X2 correlators, a subset of 
six bases is considered, as shown in Table |TT] correspond- 
ing to a choice of 3 smearings to xi and the same 3 for 
X2- For the 8x8 matrices of X1X2 correlators, a subset 
of seven bases are considered and are given in Table IIIII 
corresponding to a choice of 4 smearings to xi and the 
same 4 for X2- The correlation matrices for xi!X2 inter- 
polators contain all combinations of correlation functions 
of Xi,X2, i-e. XiXi, X1X2, X2X1 and X2X2- 

It is noted that basis operators that are linearly de- 
pendent will cause the eigenvalue analysis to fail as there 
will be a singularity in the correlation matrix. The fact 
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that our analysis (xiXi ^-nd X1X2) succeeds indicates that 
our choices of operators access an equal number of di- 
mensions in the Hilbert space. Thus it is interesting to 
examine the stability of the masses to different choices 
of bases to ascertain whether one has reliably isolated 
single eigenstates of QCD. The relevant issues are: (i) 
whether or not the operators are sufficiently far from 
coUinear that numerical errors do not prevent diagonali- 
sation of the correlation matrix and, (ii) whether or not 
the states of interest have significant overlap with the 
subspace spanned by our chosen sets of operators. Since 
our correlation matrix diagonalisation succeeded, except 
at large Euclidean times where statistical errors domi- 
nate, we conclude that our operators are sufficiently far 
from collinear. 



TABLE I: The bases of 6 x 6 correlation matrices of XiXi- 



Sweeps — >■ 


1 


3 


7 


12 


16 


26 


35 48 


Basis No. 4- 


Bases 


1 


1 


3 


7 


12 


16 


26 




2 


1 


3 


7 


12 


16 




35 


3 


1 


3 


7 




16 


26 


35 


4 


1 


3 




12 


16 


26 


48 


5 


1 




7 


12 


16 


26 


35 


6 




3 


7 


12 


16 


26 


35 



TABLE IL The bases of 6 x 6 correlation matrices of xiX2- 



Sweeps — >• 


1 


3 


7 


12 


16 


26 


35 


48 


Basis No. 4- 


Bases 


1 


1 








16 






48 


2 




3 




12 




26 






3 




3 






16 






48 


4 






7 




16 




35 




5 








12 


16 


26 






6 










16 


26 


35 





TABLE III: The bases of 8 x 8 correlation matrices of XiX2- 



Sweeps — 7> 


1 


3 


7 


12 


16 


26 


35 


48 


Basis No. 4- 


Bases 


1 


1 




7 




16 




35 




2 






7 


12 


16 


26 






3 




3 




12 




26 




48 


4 






7 


12 




26 


35 




5 






7 




16 


26 


35 




6 






7 




16 




35 


48 


7 








12 


16 


26 


35 





IV. RESULTS 
A. Variational Analysis 

We begin by considering the 6x6 correlation matrices 
of XiXi correlators. In Fig.[l] masses from the projected 
correlation functions and eigenvalues for the 3"^ basis of 
Table |T] arc presented. The basis elements for the 6x6 
matrices arc highly linearly dependent, and this means 
the analysis is more complicated than for the 4x4 ma- 
trices of Ref. [131. While the masses from the eigenval- 
ues again display larger dependency on the variational 
parameters, as observed in Rcfs. [3, [HI, masses from 
the projected correlation functions are very consistent on 
istart (= t) and At. Notably, the masses for the ground 
and first excited states are as robust as in Ref. [2J] . For 
the fifth excited state at the heaviest pion mass with the 
variational parameters (tstart, At) = (6,1),(6,2) and (7,1), 
(top left graph of Fig. [J), acceptable fits [2^ were unob- 
tainable. Nonetheless, the consistency of the calculated 
masses over the significant sets of variational parameters 
is self evident, as is how a mass can be exposed using one 
set of tstart and At. 

From a series of tstart and At, a single mass is se- 
lected for one set of tg^art and At by the selection crite- 
ria discussed in Ref. [25| , where we prefer larger value of 
istart + At [131 . In cases where a larger tgtart + At provides 
a poor signal-to-noise ratio, for example (tgtart, At)=(7,3) 
(top left figure), we prefer a little lower tgtart + At value, 
for example (tstart, At)=(7, 2), and we follow this proce- 
dure for each quark mass. 

In Fig. [21 the masses for all the 6x6 correlation ma- 
trices of XiXi correlators arc presented. Straight lines 
are drawn to illustrate the invariance of the masses over 
the bases. Since the l'^' basis contains all the consecutive 
smearing sweep counts, i.e. 1,3,7,12,16,26 sweeps, it may 
not span the space as well as other choices. This basis 
exposes higher first and second excited states and this 
remains true for the 2"*^ basis where it contains consec- 
utive five lower smearing sweeps (left figure). These are 
therefore less reliable bases for extracting excited state 
masses. It is noted that masses from the third excited 
state onwards contain higher fluctuations than the lower 
states, independent of the choice of basis. A basis de- 
pendence is also evident for the third excited state (right 
figure). However, the bases from the 3"^"^ to the 6**^ sets 
are more consistent, and in particular, the lower three 
states are stable over these four bases. Therefore, the 
bases from 3'''^ to 6"^ are used to perform a systematic 
analysis to calculate the systematic errors associated with 

the choice of basis, using at, = \l nI-i Y^f=ii^'^i " -^^)^, 
where Ni,~A in this case. 

In Fig. [31 the final results for the 6x6 correlation 
matrix analysis of XiXi correlators are shown. As the 
highest excited state accommodates remaining spectral 
strength, this state may not be an eigenstate and is not 
shown in the figure. Masses are averaged over the four 
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FIG. 1: (color online). Masses of the nucleon, states, from the projected correlation functions as shown in Eq. [32] (left) 

and from the eigenvalues (right) for the 3^'' combination (1,3,7,16,26,35 sweeps) of 6 x 6 correlation matrices of xiXi correlation 
functions. The figure corresponds to pion masses of 797 MeV (top row) and 249 MeV (bottom row). Each pair of ground and 
excited state masses correspond to the diagonalization of the correlation matrix for each set of variational parameters tstart 
(shown in major tick marks) and At (shown in minor tick marks). Here, fstart = t, as shown in Eqs. (|29p and (|31|l . 
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FIG. 2: (color online). Masses of the nucleon, states, from the projected correlation functions for the 6x6 correlation 

matrices of xiXi correlators, for pion masses of 797 MeV (left) and 295 MeV (right). Numbers on the horizontal scale correspond 
to each basis of 6 x 6 matrices, for instance, 1 and 2 correspond to the bases of (1,3,7,12,16,26 sweeps) and (1,3,7,12,16,35 
sweeps) respectively, and so on. Masses are extracted according to the selection criteria described in the text and in Ref. [2^ 
from all the combinations of 6 x 6 correlation matrices as shown in Tabled Straight lines are drawn to illustrate the invariance 
of the masses over the bases. 



bases (from the 3'^ to 6**^), whereas the errors shown are 
a combination of average statistical errors (ct^) over the 
four bases and systematic errors due to the basis choice 
(cTb). As expected from the 4x4 analysis [l^l, a simi- 
lar lower lying Roper state is also revealed in this 6x6 
correlation matrix analysis, which also has a tendency to 
approach the physical state. This analysis presents six 



distinct energy states, where the ground and first excited 
states are extremely robust, emphasizin g th e utility of 
the analysis technique presented in Ref. |24| . From this 
figure it is also evident that the second excited state re- 
mains high and an unnatural ordering of the error bars 
is manifest. Hence there are concerns as to whether this 
state is a true energy state. 
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FIG. 3: (color online). Masses of the nucleon, N2 states, for the ground (g.s) and the excited (e.s) states from the 6x6 
correlation matrices of XiXi correlators over four bases (from the 3"''* to the 6* of Fig. [2]). The errors shown in the figure are a 
combination of average statistical errors over these four bases and systematic errors due to basis choices. Errors are combined 
in quadrature. The black filled symbols are the experimental values of the ground and the Roper states obtained from Ref. [s^. 
Masses are given in Table [TVl 



TABLE IV: Masses of the nucleon, states, averaged over the four bases (from the a"'* to 6"'). The errors shown here are a 

combination of average statistical errors over four bases and systematic errors for the choice of basis, combined in quadrature. 





aMg.s 


aMist (,.s (Roper) 




aMyd 




0.5141(19) 


1.0414(68) 


1.451(40) 


1.784(93) 


2.447(60) 


2.813(76) 


0.4705(20) 


0.9933(78) 


1.409(41) 


1.892(98) 


2.459(49) 


2.856(78) 


0.4134(22) 


0.9286(80) 


1.357(43) 


1.875(54) 


2.451(54) 


2.838(90) 


0.3490(24) 


0.8588(88) 


1.305(45) 


1.857(61) 


2.436(41) 


2.849(59) 


0.2776(24) 


0.781(10) 


1.250(53) 


1.886(65) 


2.427(39) 


2.845(54) 


0.2452(24) 


0.752(11) 


1.229(60) 


1.903(38) 


2.423(42) 


2.846(58) 


0.2110(27) 


0.722(14) 


1.204(75) 


1.894(42) 


2.417(45) 


2.850(72) 


0.1905(31) 


0.711(12) 


1.180(89) 


1.892(46) 


2.414(47) 


2.840(74) 


0.1607(35) 


0.682(14) 


1.142(11) 


1.891(51) 


2.410(55) 


2.837(72) 



In Rcf. [36[ addressing negative parity nuclcons, a 
whole new set of nearby states is observed upon intro- 
ducing X2- Thus, it is important to explore the role 
of X2 in the positive-parity sector. In Fig. [4l masses 
from the projected correlation functions and eigenvalues 
for the 6x6 correlation matrix analysis of X1X2 corre- 
lators are presented. The diagonalization is only suc- 
cessful for the heavier four quarks, which again proves 
that a variational analysis with X2 interpolator is always 
a challenge [2^. This interpolator vanishes in the non- 
relativistic limit 0, [s^l and is renowned for providing 
noisier correlation functions [ll|, [Hi especially for the 
lighter quark masses, where one otherwise might expect 
the relativistic nature of the quarks to be of benefit. The 
Euclidean-time correlation functions of the X2 interpola- 
tor decay more rapidly, and the signal-to-noise ratio also 
deteriorates more quickly with time. In Refs. [13, [25j . 
a rise from below of excited state masses from the pro- 



jected correlation functions was observed for the linearly 
dependent operators. This is also manifest in the 6x6 
analysis of xiXi, shown in Fig. [H but for the 6x6 anal- 
ysis of X1X2 correlators, these effects are reduced. This 
is expected because the basis elements of these bases are 
likely to be more orthogonal, improving the isolation of 
states. 

In Fig. m projected masses from all the 6x6 correla- 
tion matrices of X1X2 ars presented. Again the ground 
state is robust in this analysis. As the 6"^ combination 
contains consecutively higher number of smearing sweeps 
(16,26,35), this basis is not as reliable as the other sets. 
The first excited state extracted with this basis sits some- 
what lower than the other five bases. A careful analysis 
of Figs. [5] and [2] reveals that in the vicinity of the second 
excited state from the 6x6 analysis of XiXi correlators 
(Fig. [2]), three excited states appear in the 6x6 analysis 
of X1X2 (Fig-Hi)- ^ completely new excited state, the sec- 
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FIG. 4: (color online). Masses of the nucleon, A'' 2 states, from the projected correlation functions as shown in Eq. [32] (left) 
and from the eigenvalues (right) for the 3'^'^ combination (3,16,48 sweeps) of 6 x 6 correlation matrices of X1X2 correlation 
functions. The figure corresponds to the pion mass of 797 MeV. Each pair of ground and excited states mass correspond to 
the diagonalization of the correlation matrix for each set of variational parameters tstart (shown in major tick marks) and At 
(shown in minor tick marks). 
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FIG. 5: (color online). Masses of the nucleon, states, from the projected correlation functions for the 6x6 correlation 

matrices of correlators, for pion masses of 797 MeV (left) and 541 MeV (right). Numbers on the horizontal scale 

correspond to each basis of 6 x 6 matrices, for instance, 1 and 2 correspond to the bases of (1,16,48 sweeps) and (3,12,26 
sweeps) respectively, and so on. Masses are extracted according to the selection criteria described in the text and in Ref. [2^| 
from all the combinations of 6 x 6 matrices as shown in Table |IT1 Straight lines are drawn to illustrate the invariance of the 
masses over the bases. 



ond excited state, is extracted in the X1X2 analysis. This 
second excited state is robust for all bases (left figure). 
The signal-to-noise ratio for this state deteriorates more 
rapidly for the light quark (right figure). Interestingly, 
the fifth excited state extracted from the X1X2 matrices, 
sits at a significantly lower energy state than those from 
the XiXi matrices, refiecting the fact that two new low- 
lying states have been revealed. The third and fourth 
excited states are nearly degenerate in this analysis. As 
the results for the lower three states are consistent over 
five basis choices (from 1**' to 5*^) of the 6x6 matrices, 
we perform a systematic analysis over these five bases, as 
discussed previously. 

In Fig. [6l a summary of the 6x6 correlation matrix 
analysis of X1X2 correlators is presented. Masses are av- 
eraged over the five bases, while errors are a combination 
of average statistical and systematic errors due to ba- 



sis choices. Apart from the ground and the first excited 
states, this analysis presents a distinct second excited 
state for the three heavier quark masses. The signal for 
the light quark mass for this state is poor. This analysis 
reveals new nearly degenerate third and fourth excited 
states, that the 6x6 bases of xiXi correlation functions 
arc unable to resolve. 

In Fig. [71 masses from the projected correlation func- 
tions and eigenvalues for the 8x8 correlation matrix 
of X1X2 correlators are presented. Similar to the 6x6 
analysis of XiX2j this analysis is also successful for the 
heavier four quark masses. The enhanced dimension of 
the matrices means the numerical diagonalization is less 
stable, and so the variational analysis is only success- 
ful for a few sets of variational parameters, tstart and At. 
However, in Fig. [71 there is sufficient consistency between 
the masses from the projected correlation functions for 
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FIG. 6: (color online). Masses of the nucleon, N states, for the ground (g.s) and the excited (e.s) states from the 6x6 
correlation matrices of X1X2 correlators over five bases (from the 1"' to S"" of Fig. The errors shown in the figure are a 
combination of average statistical errors over these five bases and systematic errors due to basis choices. Errors are combined 
in quadrature. Masses are given in Table [V] 



TABLE V: Masses of the nucleon, A'' 2 states, averaged over the five bases (from the 1=* to 5*''). The errors shown here are a 
combination of average statistical errors over these bases and systematic errors for the choice of basis, combined in quadrature. 





aMg.s 


aM^st (Roper) 








0.5141(19) 


1.0419(65) 


1.526(28) 


1.669(37) 


1.960(64) 


1.85(12) 


0.4705(20) 


0.9927(69) 


1.474(39) 


1.638(45) 


1.920(72) 


1.84(12) 


0.4134(22) 


0.9295(75) 


1.432(33) 


1.617(66) 


1.875(79) 


1.82(12) 


0.3490(24) 


0.8612(84) 


1.377(39) 


1.63(11) 


1.81(12) 


1.74(17) 
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FIG. 7: (color online). Masses of the nucleon, states, from the projected correlation functions as shown in Eg. 1321 (left) 

and from the eigenvalues (right) for the 3'^'^ combination (3,16,48 sweeps) of 8 x 8 correlation matrices of X1X2 correlation 
functions. The figure corresponds to the pion mass of 797 MeV. Each pair of ground and excited states mass correspond to 
the diagonalization of the correlation matrix for each set of variational parameters tstart (shown in major tick marks) and At 
(shown in minor tick marks). 



each set of tstart and At. As before, the 8x8 analysis 
gives a projected mass which is highly independent of the 
variational parameters. Here we have selected istart = 6, 
At = 2. 



In Fig. in masses from the projected correlation func- 
tions for all combinations of 8 x 8 correlation matrices 
are shown. All the bases reveal a very consistent ground 
state mass. A systematic basis dependency of the excited 
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FIG. 8: (color online). Masses of the nucleon, states, from the projected correlation functions for the 8x8 correlation 

matrices of X1X2 correlators, for pion masses of 797 MeV (left) and 541 MeV (right). Numbers on the horizontal scale correspond 
to each basis of 8 x 8 matrices, for instance, 1 and 2 correspond to the bases of (1,7,16,35 sweeps) and (7,12,lfL26 sweeps) 
respectively, and so on. Masses are extracted according to the selection criteria described in the text and in Ref. 25] from all 
the combinations of 8 x 8 matrices as shown in Table IIIII Straight lines are drawn to illustrate the invariance of the masses 
over the bases. 
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FIG. 9: (color online). Masses of the nucleon, 
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states, for the ground (g.s) and the excited (e.s) states from the 8 x 



correlation matrices of X1X2 correlators over four bases (from the 3"^"* to fi"^ of Fig. [8]). The errors shown in the figure are a 
combination of average statistical errors over these four bases and systematic errors due to basis choices, combined in quadrature. 
Masses are given in Table IVll 



TABLE VI: Masses of the nucleon, states, averaged over four bases (from the 3'''* to 6*''). The errors shown here are 

a combination of average statistical errors over these four bases and systematic errors for the choice of basis, combined in 
quadrature. 





aMg.s 


aMi^t ,,,(R) 


aMjnd 






aMgth c.s 




0.5141(19) 
0.4705(20) 
0.4134(22) 
0.3490(24) 


1.0412(66) 
0.9912(71) 
0.9277(77) 
0.8588(86) 


1.470(33) 
1.427(35) 
1.372(38) 
1.311(42) 


1.669(33) 
1.642(38) 
1.619(52) 
1.623(92) 


1.874(78) 
1.825(79) 
1.758(77) 
1.779(80) 


1.843(90) 
1.888(80) 
1.896(48) 
1.895(57) 


2.300(67) 
2.309(54) 
2.304(45) 
2.269(52) 


2.478(58) 
2.475(46) 
2.462(47) 
2.475(41) 
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states (from third excited state onwards) is noticed for 
the I''* and 2"'^ combinations. As the l'^' combination 
consists of 1 and 7 smearing sweeps, it provides a higher 
mass from the third excited state onwards, in comparison 
with the other bases. The 2"'^ basis contains consecutive 
smearing sweep counts, starting with 7 sweeps, and end- 
ing with a moderate sweep count of 26. This provides 
less diversity in this basis for the higher excited states 
and also provides a higher mass from the third excited 
state onwards. Therefore, the 1*^* and 2"'^ bases are not 
as reliable as other sets. The fifth and sixth excited states 
sit a little high for the 7"^ basis which has a narrow clus- 
ter of smearing sweeps and is therefore likely to be less 
reliable in spanning the space. Bases from the 3'^'^ to 6*^ 
sets have more smearing diversity being well spread over 
the range of smearing sweep counts starting from 3 and 7 
sweeps. It is also evident from Fig. [8] that the S'*^ to 6*"^ 
combinations of smearing sweeps are successful in pro- 
viding highly consistent lower energy states. As before, 
we will therefore perform a systematic analysis over these 
four bases. It is interesting to note that in the vicinity 
of the fifth excited state from the 6x6 analysis of X1X2, 
the 8x8 correlation matrix presents three excited states. 

In Fig. [21 a summary of the 8x8 correlation matrix 
results are presented. As in Fig. [51 a distinct second ex- 
cited state is also obtained in this 8x8 analysis, for the 
two heavier quark masses. However, increased separation 
of the third and fourth excited states appears here with 
the enlarged basis. A careful examination reveals that 
the first excited state extracted in this analysis, is a little 
lower than in the 6x6 analysis of X1X2 (see Fig. [6]). In 
accord with the principle of the variational method, our 
analysis signifies the importance of using larger correla- 
tion matrices to reliably isolate the higher energy eigen- 
states of QCD. 



B. Splitting of Excited States 

In Fig. [TUl we show how the excited energy states de- 
pend on the choice of correlation functions and the di- 
mension of the correlation matrix. The ground state is 
clearly independent of the dimension of the matrices and 
choice of interpolating fields. The excited state extracted 
from the 2x2 correlation matrix splits into two excited 
states with the 3x3 matrix. The central value of the 
first excited state, extracted with the 3x3 matrix, is a 
little lower than that of 2 x 2 case. Interestingly, the sec- 
ond excited state from the 3x3 matrix is in very good 
agreement with the excited state number of the 2x2 
analysis. 

The 4x4 analysis of XiXi correlation functions reveals 
a lower first excited state (Roper state) [l^l and two other 
heavier excited states. Therefore, the 4x4 basis of xixi 
correlation functions is able to resolve and isolate the su- 
perposition of eigenstates in the 3x3 matrix. The 6x6 
and 4x4 correlation matrices of xiXi correlators provide 
very consistent results for the lower four energy states. 



> 

S3 




Various matrices 



N^'*' states. 



FIG. 10: (color online). Masses of the nucleon 
for the various dimensions of correlation matrices, as labelled 
on the upper horizontal axis. The 2x2 and 3x3 results for 
the smeared-smeared correlation functions with 26 sweeps of 
smearing are taken from Ref. [2^, while the 4x4 results are 
from Ref. [Ij. The 6x6 and 8x8 results correspond to 
the analyses presented in this paper. Dotted lines are drawn 
to aid in illustrating the consistency of the results. Figure 
corresponds to the pion mass of 797 MeV. 



and the 6x6 analysis is able to extract two new higher 
energy states. The larger error bar of the third (high- 
est) excited state extracted with the 4x4 correlation 
matrix compared to the 6x6 result of xiXi is a manifes- 
tation of accommodating remaining spectral strength in 
this (third excited) state in the 4x4 analysis and hence 
is unreliable [23 |. 

The 6x6 analysis of x 1X2 correlation functions extracts 
different excited states than those of the 6x6 analysis of 
XiXi alone. The first excited state extracted with this ba- 
sis sits high in comparison with the Roper state. Thus, a 
basis of four different smearings with XiXi is essential to 
isolating the single eigcnstate associated with the Roper 
resonance. Another interesting feature is the splitting of 
the "second excited state" of the XiXi analysis into three 
nearby states with the new X2 spin-flavour combinations. 

The final very interesting outcomes are revealed from 
the 8x8 correlation matrix analysis of xiX2- This anal- 
ysis provides a first excited state which is in excellent 
agreement with the Roper state revealed by XiXi corre- 
lators alone. This present investigation reveals that X2 
plays a subtle role for the Roper. We note however that 
a new nearby second excited state is revealed in the X1X2 
analysis and a future high-statistics analysis may reveal 
it is essential to resolve this state correctly to obtain the 
first excited state accurately. This basis also provides sec- 
ond, third and fourth excited states consistent with the 
6x6 analysis of XiX2- However, the fifth excited state 
coming from this analysis is a little lower than that of 
the 6x6 matrix of XiX2- There the fifth state accommo- 
dates all remaining spectral strength. This state wasn't 
identified in the XiXi analysis. The sixth excited state 
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FIG. 11: (color online). Masses of the nucleon 
from the 8x8 correlation matrices of X1X2, and 4x4 correla- 
tion matrices of XiXi ^ ^ function of the squared pion mass, 
for the ground (g.s), first through sixth excited states (e.s). 
The experimental values are taken from Ref. [s^, where the 
reliability of the ground and Roper states are signified by four 
stars (****), Pu (1710 MeV) state by three stars (***) and 
Pii (2100 MeV) state by one star (*). 



provided by this analysis is consistent with the third ex- 
cited state extracted with the 6x6 analysis of XiXi- It is 
worth noting that the basis of XiXi correlation functions 
is insufficient to isolate the second excited state of the 
nucleon. 

While this analysis is able to provide some evidence 
that the first five energy-states are reasonably robust 
through comparisons of 6 x 6 and 8x8, only a more com- 
prehensive higher dimension correlation matrix analysis 
can assess the reliability of the sixth and seventh states. 

In Fig. Ill) a summary of results for the positive p arity 
excited states of the nucleon is presented. In Ref. j35| . 
the quality of the Pn (1710) MeV state is characterized 
by three stars (***) and Pn (2100 MeV) state by one star 
(*). These resonance states decay through the N — >■ Ntttt 
channel with positive parity. Therefore, looking at both 
the S-wave N + tt + Tr and P-wave N + tt decay channels 
provides information related to these resonances. The 
second excited state is very close to the threshold N + n 
state but has a different slope. It is also evident that this 
state is significantly different and lower than the iV-t-7r-f tt 
state for the heavier three quark masses providing evi- 
dence that this state is best described as a single particle 
state. It is also noted that the third and fourth excited 
states have a mass dependence different from the mul- 
tiparticle states. This suggests that these two energy 
states revealed here are more likely to be single particle 
resonance states. Simulations on larger lattice volumes 
and at lighter quark masses are necessary to resolve this 
issue. 

Fig- HH provides a comparison of the nucleon spectrum 
revealed in this comprehensive analysis with state-of-the- 



art results from Refs. [l7[ and [38| in quen ched QCD. 
The results from the BGR collaboration [3^ follow from 
an analysis similar to ours where fermion source smear- 
ing is key to providing a basis which enables the isola- 
tion of excited states. However, they use only two levels 
of source smearing and use different smearings for each 
quark fiavour within a given interpolator. The fact that 
their Roper state sits high relative to our result is con- 
sistent with our observation that a minimum of four dif- 
ferent smearing levels is required to obtain a stable low- 
lying state. Whereas our 6x6 analysis with six different 
smearing levels did not change the mass of the Roper 
state from our 4x4 analysis, the consideration of only 
three smearing levels did lead to a higher mass Roper 
state. 

Results from the LHP collaboration [l3] are illustrated 
at TOtt = 490 MeV. Here we have reported all states that 
may have spin-1/2 quantum numbers. The Roper state 
in their analysis sits higher than both the present analysis 
and the results of Ref. [Ill ■ their scale determination 
provides a low nucleon mass, the excited state masses will 
become higher again if one sets the scale to reporoduce 
the nucleon mass of this analysis or that of Ref. [11] . It 
suggests that their consideration of a single level of source 
smearing may be insufficient to provide a basis suitable 
to resolve the Roper state. Our expectation is that the 
introduction of a variety of source smearings into their 
basis operators may be beneficial. It would be interesting 
to expand their basis to include several levels of source 
smearing to explore the extent to which the excited state 
masses change. 



V. CONCLUSIONS 

In this paper, we have presented a comprehensive cor- 
relation matrix analysis of the positive parity spin-i exci- 
tations of the nucleon. We have considered large dimen- 
sions of correlation matrices, built from smeared-smeared 
correlation functions of xixi and X1X2 correlators. The 
results of this paper, for the ground and Roper states of 
the nucleon, are very consistent with the discovery of the 
Roper in Ref. [2^ from a 4 x 4 correlation matrix of XiXi 
alone. Thus, the most important conclusion is that the 
earlier first result for the Roper [23| is robust under the 
comprehensive analysis presented herein. 

This paper signifies the importance of using a variety 
of smearings in making an operator basis and the impor- 
tance of comparing the results across a range of choices 
of basis for the correlation matrix. In particular, it em- 
phasizes the use of at least four different smearing levels 
to isolate the elusive Roper state. 

In addition, we have reported new second, third, 
fourth, fifth and sixth excited-states of the nucleon. 
While both the XiXi and X1X2 correlation matrices pro- 
vide excellent agreement for the ground and the Roper 
states, only the X1X2 analysis provides access to a com- 
pletely new (second) excited eigenstate. All the correla- 
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FIG. 12: (color online). Masses of low-lying nucleon states from this work (legend at top right), from Ref. [l7| (legend at 
top left) and Ref. [s^ (legend at bottom right). The results from Ref. [l3| correspond to = 490 MeV, where the nearly 
degenerate fourth energy state has been shifted a little to the right for clarity. 



tion matrix analyses provide a very consistent third ex- 
cited state, see Fig. [TOl Simulations on larger lattice vol- 
umes at lighter quark masses with higher statistics will be 
interesting to investigate the propagation of these states 
towards light quark masses. 

Another interesting result of this paper is delineating 
the nature of the splitting of the excited states of the 
nucleon, as seen in Fig. [TOl for the variational method. 
While the ground state is independent of any basis choice, 
and the dimension of the matrix used, extracting a partic- 
ular excited state can depend on the interpolators mak- 
ing the correlation functions. For example, a variational 
analysis of XiXi correlation functions doesn't provide the 
second excited state, indicating that the xi interpolator 
either doesn't couple or only has a very small coupling 
to this state. We emphasize the importance of using a 
large correlation matrix for a reliable extraction of ex- 
cited states energies. 

Because of the difhculties in ascertaining that a given 
set of basis interpolators is sufficient to isolate all the 
eigenstates of QCD in the mass range of interest, we be- 



lieve that it is essential to resolve the nucleon spectrum 
using a wide variety of approaches. For example, it is 
now apparent that the consideration of a single level of 
fermion source smearing is unlikely to be sufficient to 
resolve the low-lying Roper state. Future investigations 
should explore the role of fermion source smearing in pro- 
viding a suite of basis operators that span the space in 
an effective manner. 

To test the robustness of the spectrum revealed in this 
analysis it will be important for future studies to carefully 
examine a wide variety of different interpolators in order 
to test whether or not other low-lying states will appear 
in the spectrum. For example, we anticipate that five- 
quark operators will be essential in accessing the relevant 
low-lying meson-baryon scattering states. 
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